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ABSTRACT 

We present the rest-frame ultraviolet through near infrared spectral energy distribution for an 
interacting Lyman break galaxy at a redshift z = 4.42, the highest redshift merging system known 
with clearly resolved tidal features. The two objects in this system - HDF-G4 and its previously 
unidentified companion - are both -B435 band dropouts, have similar Vgoe — *775 and 1775 — zgso 
colors, and are separated by 1", which aX z — 4.42 corresponds to 7 kpc projected nuclear separation; 
all indicative of an interacting system. Fits to stellar population models indicate a stellar mass of 
Mi, — 2.6 X lO^^M©, age of r^, = 720 My, and exponential star formation history with an e-folding 
time To — 440 My. Using these derived stellar populations as constraints, we model the HDF-G4 
system using hydrodynamical simulations, and find that it will likely evolve into a quasar by z ~ 3.5, 
and a quiesent, compact spheroid by z ^ 2.5 similar to those observed at z <; 2. And, the existence 
of such an object supports galaxy formation models in which major mergers drive the high redshift 
buildup of spheroids and black holes. 

Subject headings: cosmology: observations - galaxies: evolution ~ galaxies: high-redshift - galaxies: 
stellar content - galaxies: interacting - infrared: galaxies 



1. INTRODUCTION 



Recently, there have been great advances in un- 
derstanding the nature and evolution of high-redshift 
alaxies. T hroug h either col o r-sele ct ion criteria 
Steidel et al.] 120031 iFranx et all 12 003': 'Dadd i et all 



20041) ■ Lyman-g emission (iRhoads 
Malhotra fc Rhoads ' 2002; ' Aiiki et all 



Malhotra 120011 
L2003; H u et all 



2004l:lTaniguchi et al... 2005.) or blank-field submillimeter 



surveys combined with radio observations (iSmail et al 



■veyi 

2OOOI: iBa rger et al.ll2000t llvison e t al."2002': 'Borv s et al 
2003; C hapman et al.l l2005t ICoppin ct al. 200£j) ob- 
servers have compiled large catalogs of z 2 galaxies. 
Subsequent photometric and spectroscopic foUowup has 
been successful at constraining the physical properties 
of these galaxy populations, including the evolution of 
the star forma tion rate (SFR) and stellar mass density 
out to z ~ 3 fShaplcv ct al. 2003; Barmbv et all 120041 : 
IShaplev et all feoOS: Lai et al. 2007). 

One method for selecting high-re dshift galaxies - the 
Lyman break dropout techinque (| Steidel fc HamiltonI 
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fl99l - uses the 91 2 A break to select z ~ 3 galaxies. 
Spectroscopic foUowup has pro vided subst antial samples 
of confirmed z ~ 3-4 galaxies (|Steidel et a l. 1999, 2003). 
These Lyman break galaxies (LBGs) are known to have 
high star formation rates relative to local galaxies of the 
same stellar ma ss, as revealed by their high rest frame 
UV luminosity (jSteidel et al.ll2003D . Until recently, stel- 
lar masses for these objects were estimated using ground- 
based optical photometry. At these wavelengths in the 
rest-frame, the luminosity is dominated by recent star 
formation, rather than the accumulated mass from older 
stellar populati ons. Studies usin g the Infrared Array 
Camera (IRAC: iFazio et al1l2004[ ) on the Spitzer Space 
Telescope, which provides photometry out to 8/xm in the 
observed frame, have used the rest-frame K band lumi- 
nosity to place more robust constrai nts on the bulk of 
the s tellar content of LBGs at z ~ 3 (jRigopoulou et al.) 
[20061) . 

At the same time, r ecent theoretical model- 
ing (iHopkins et aTl 2006a, b, l2007al lbl) and simulations 
(iCox et al.ll2006F have suggested a link between merg- 
ing galaxy populations, quasars, and present day el- 
lipticals through the self-regulated growth of super- 
mass ive black holes f SMBHs) in gas-rich major merg- 
ers (jPi Matteo et all [2005.1 T he pr e sence of large 
amounts of dus t (ISawicki fc Yed 119981 ICalzettil 120011 : 
iTakeuchi fc Ishiill2004r ) that implies a high SFR, m com- 
bination with clustering arguments (jAdelberger et aTl 
|2005( ) , and a me rging fraction of ^ 10 — 25% at high red- 
shift (jLotz et a l. 2006) have suggested empirically that 
high-redshift LBGs ma y be the progenitors of those same 
present-day ellipticals (jPettini et al.lll998[ ). As a result, 
major mergers involving LBGs at high redshift may pro- 
vide observational evidence of a link between these pop- 
ulations. 

In this work we present optical through infrared ob- 
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servations of HDF-G4, an LBG with a spectroscopi- 
cally confirmed redshift of z — 4.42 in the ffubble Deep 
Field North (HDFN) with a(J2000) = f2: 37:20.57 and 
( 5(J200 0) =+62:ff:6.08 from the catalog of ISteidel etlHI 
()1999t ). This object, and its previously unidentified 
companion, constitute the most distant interacting sys- 
tem observed to date for which the resolved morphology 
clearly suggests a merger. In § [2l we present multiwave- 
length data for the HDF-G4 system. From these data, in 
§ [21 we use stellar population synthesis models to derive 
stellar masses and star formation rates (SFRs) for each 
component. In § HI we use the observed stellar popula- 
tions and morpholgy to constrain a model of the HDF- 
G4 system using hydrodynamical simulations. Finally, 
in § [5l we discuss the implications of the predictions of 
our model for the future evolution of the IIDF-G4 sys- 
tem. Our analysis suggests that HDF-G4 will be a high- 
redshift passive spheroid by z ~ 2.5, similar to thos e ob- 
served bv lLabbe etHl (|2005f ) and lZirm etlll (|2007f ). For 
the proceeding analysis, we assume the concordance cos- 
mological model; a flat ACDM cosmology with JIa = 0.7 
and h = 0.7. All magnitudes presented are in the AB 
system. 

2. MULTI- WAVELENGTH PHOTOMETRY OF HDF-G4 

The Great Observatories Origins Deep Survey 
(GOODS), which includes the HDFN, is the deepest 
multi-wavelength survey, detecting gala xies at extremely 
high redshifts iii optical and IR bands (|Mobasher et al.1 
120051 lYan et all [200l iLai et all 120071 ). and potentially 
interacting sy stems at very h i gh re dshift at z ^ 4 
(|Rhoads et all l2005t lYan et al.ll2005l ). The wealth of 
multi-wavelength data available for GOODS fields, from 
the X-Ray to radio, makes it ideal for studying the prop- 
erties of such interesting systems. 

We visually inspected the available optical imaging 
data, and find that HDF-G4 appears to be one ob- 
ject in the ground-based Subaru images (jCapak et al.1 
I2004D. but is resolved in to two objects in the ACS mosaic 
(|Giavalisco et al.ll2004D . Aperture photometry obtained 
with Sextractor, a publ icly available source dete ction 
and photometry package (jBertin fc Arnouti Il996f ). for 
both the primary (CI) and secondary (C2) components 
were taken from the public catalog for all four ACS mo- 
saics: F435W (B435), F606W (T^soe), F775W (^775), and 
F850LP (zsso). C2 is 1 magnitude fainter in the ACS 
1775 band, and ^ 1" away from CI. Furthermore, the 
position angle and center of the slit-mask indicate that 
the observed optical spectrum is of CI only. 

There are three pieces of evidence suggesting that 
CI and C2 are physically associated with each other. 
The first is that both objects are i?435 dropouts with 
similar observed optical colors, indicating that they 
are at comparable redshift. The second is their rel- 
ative angular separation, which a,t z = 4.42 corre- 
sponds to a projected internuclear distance of 7 kpc. 
This is simil ar to local interacting systems such as the 
"Antennae" (IWhitmore fc Schweizeil HOOl . II ZW 96 
(jColdader et al.lll997D. or ultraluminous infrared galax - 
ies (ULIRGs HRigopoulou et al.lll999t ISurace et all2000D . 
The third, and most compelling piece of evidence is a ma- 
terial bridge connecting both components. This feature 
is low surface brightness - typically no more than 2a 
above the sky noise per pixel - but highly statistically 




Fig. 1. — Stacked optical imaging data centered on HDF-G4 
for all four ACS filters: -B435, Vaoa, 1775, and zg^Q. Red contours 
correspond to 2-5 standard deviations above the noise. The green 
circle has a 1" radius, which at this redshift corresponds to ~ 7 kpc 
in our assumed cosmology. Note the material exchanged between 
the two sources, contained within a 2— cr contour. This is, despite 
its low per— pixel surface brightness, a highly statistically significant 
feature that indicates an interacting system. 
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Fig. 2. — Infrared imaging data, with negative contrast, centered 
on HDF-G4, in all four IRAC channels: (left to right) 3.6, 4.5, 5.8, 
and 8.0 fj,m. The red circles have a 1" radius. HE)F-G4 is clearly 
detected at 3.6/im and A.bfim , and marginally detected at 5.8/im 
and 8.0^m . 

significant; a Monte-Carlo analysis indicates that it is 
detected at lOtr in the stacked optical image (see Fig- 
ure [T|). It furthermore has similar colors to CI and C2, 
and is not detected in ^435 (see Table [1]). 

This system is also marginally detected in the CFHT 
WIRCam J and K band images (Lin 2006, Simard 2006, 
private communication). The centroid of the K-band 
counterpart is 0.35" offset from CI, but 0.8" from C2. 
Given the 0.9" seeing in the J and K images, we argue 
the detected J- and K-band flux is dominated by CI. 

The GOODS IRAC image is very deep, and thus both 
CI and C2 should be detected at least at 3.6 and 4.5/xm. 
Indeed, we find significant detections at 3.6 and 4.5/im, 
and marginal detections at 5.8 and 8.0/im (see Figure[2]). 
We obtain aperture photometry (see Table [1]) for HDF- 
G4, after subtracting out bright neighbors to minimize 
contamination using StarFinder, a publically availabl e 
PSF-fitting photometry package (jDiolaiti et al.l [lOOOf) . 
Errors were estimated using a Monte-Carlo analysis. 
IRAC, however, has 2.1" resolution, and thus cannot re- 
solve CI and C2. We use the zgso flux ratio to estimate 
the C2 contribution to the total IRAC flux of this system. 
C2 has a somewhat bluer color with Vgoe "-^sso = 0.64, as 
compared to CI with Veoe — ^sso = 1-09. Th is is very sim- 
ilar to a merging system at z=3.01 found bv lHuang et al.l 
(|2006( l. in which the IR flux densities are dominated by 
the red luminous component. We argue that C2 should 
also have bluer zgso — [3.6] color than CI, and use the 
Z850 flux ratio to estimate an upper limit of 10% on the 
contribution of C2 to the unresolved IRAC flux densities. 

The HDFN was additionally covered b y the Very Larg e 
Array at 1.4 GHz to a depth of 40 /iJy i)Richardsil2000^ . 
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TABLE 1 

Photometry of Different Components of the HDF-G4 System 





B435 




4775 


^850 


J 




3.6/im 


4.5/im 


5.8/im 


8.0/im 




[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 






[MJy] 




Cl 


> 27.2 


25.87 ± 0.02 


24.72 ± 0.02 


24.77 ± 0.02 


23.65 ± 0.59 


23.59 ± 0.18 


1.34 ± 0.22 


0.92 ± 0.34 


0.68 ± 3.47 


0.79 ± 0.94 


C2 


> 27.2 


26.43 ± 0.03 


25.67 ± 0.03 


25.80 ± 0.04 














Bridge^ 


> 27.74 


27.92 ± 0.19 


26.78 ± 0.12 


26.72 ± 0.15 
















0.1 ; 

Rest-Frame X (/im) 



Fig. 3. — The SED of HDF-G4 (open circles) and its companion 
(open triangles) from the rest-frame UV to the NIR, along with 
the best-fit Bruzual & Chariot (2003) stellar population synthesis 
models for an exponential (red; ESF) and continuous (green; CSF) 
star formation history. IRAC data for 5.8/xm and S.O/^m are shown 
for completeness, but are excluded from the fitting because of their 
low significance. 

MIPS to a depth of ~ 70/iJy at 24/im , and the Chandra 
X-Ray Telescope to depths of 2.5 x 10"^^ and 1.4 x lO^^^ 
erg cni-2 in the sof t (0.5-2.0 keV) and ha rd (2-10 keV) 
bands respectively (|Alexander et al.l I2003D . There was, 
however, no detection of either HDF-G4 or its compan- 
ion in the publicly released source lists or imaging data. 
This is all consistent - and in particular the lack of 24/im 
or 1.4 GHz emission - with the high redshift of HDF-G4. 

We present the spectral energy distribution (SED) of 
HDF-G4 (open circles) and and its companion (open tri- 
angles), along with best fit models (see §[31) in Figure O 
All the available photometry for Cl, C2, and the bridge 
are presented in Table [TJ 

3. ANALYSIS OF THE STELLAR POPULATIONS 

Our observations of each component of the HDF-G4 
system can be used to constrain the underlying stellar 
population at the time of our observations. To estimate 
the age and mass of the stellar population of the pri- 
mary component of HDF-G4 (Cl), we fit its SED to a 
gri d of population synthesis models (Bruzual fc CharlotI 
|2003), assuming a Salpctcr (1955) IMF, with three dif- 
ferent parameterized star formation histories: a single 
stellar population (SSP) instantaneous burst, continuous 
star formation (CSF), and an exponential decay r-model 
(ESF). We consider solar and sub-solar { Zp,/2QQ) mod- 
els, in addition to includin g the effects o f [Calzetti et al.l 
([2000) dust extinction and lMadaul ([1^995) Lyman-a for- 
est absorption. Our fits exclude observations in IRAC 
Channels 3 and 4, which are poorly constrained by the 
measurement. 

Results for the best -fit CSF and ESF models are shown 
in Figure [21 the SSP models were excluded because they 
were both a poor fit to the data. Requiring the age 
of the stellar population to be less than a Hubble 



time at z = 4.42 further excludes the sub-solar CSF and 
ESF models. The best fit overall, with xVd-o.f. = 1.46 
was the solar metallicity ESF model, with e-folding time 
To = 440 My, stellar age r^, — 720 Myr, stellar mass 
A'h = 2.6 X 101°Mq, and dust extinction E{B -V) ^Q. 
The best-fit CSF model, with xVd.o.f. = 2.24, had 
= 1020 Myr, NU = 2.6 x IO^^Mq, E{B -V) ^ 0.05, 
and star formation rate SFR = 32 M© yr~^ . These re- 
sults are so mewhat s e nsitiv e to the choice of IMF; e.g., 
assuming a iChabriei) (|2003D distribution will reduce the 
inferred stellar masses by ~ 30—40%. These inferred stel- 
lar masses and SFRs are broadly consistent with stud- 
ies applyi ng a similar analy s is to large populati ons of 
LBGs (Papovich etld] [200]] : IShaplev et al.l l200lD . and 
with 8/im selecte d LBGs with similar infr ared-to-optical 
colors at z ~ 3 ("Rigop oulou et al.]|2006[ ). We note that 
our best fit stellar populations indicate minimal dust ex- 
tinction, which is unusual for LBGs but is consistent with 
the least obscured object s at z 2 — 3 (jPapovich et all 
[200lt IShaplev et al.l[200ll) . 

We do not have sufficient spectral coverage to fit mod- 
els to the photometry for C2 or the bridge. However, 
at z = 4.42, the observed optical images cover the rest- 
frame ultraviolet continuum from 800-2310A. The inte- 
grated spectrum in this waveband is dominated by young 
stars, an d there fore scales linea rly with the star forma- 
tion rate (jKenni cutt 1998; Mada u et a"Llll998f) . Using the 
calibrations of [Kennicuttl (|199 8) for the intrinsic luminos- 
ity at 1500 A, which assume CSF and a lSalpetej ([l955i ) 
IMF, we find that our 1775 (rest-frame Ae// « 1421 A) 
photometry implies SFR « 9 Mq yr~^ . Similarly, the 
bridge traces significant off-nuclear star formation, with 
SF R 3 M(7i yr ~^ as inferred from the 1775 magnitude. 
If a IScalol (|1986f ) IMF is used, these estimates increase 
by a factor of ~ 2. 

The lack of a detection of HDF~G4 in deep Chandra 
observations can constrain the relative contributions of 
an active galactic nucleus (AGN) and stellar components 
of the bolor netric luminosi t y in Cl . From the AGN SED 
template of lHopkins eTall (I2007dl ). we find that the lim- 
iting fiux of the hard band (2 — 8 keV) puts a strict 
upper limit, in the absence of significant scattering or 
absorption by dust, of Lagn ^ 2 x lO^^L© on the bolo- 
metric AGN luminosity. We believe that neglecting the 
effects of dust on the emitted X-rays is a good approx- 
imation for two reasons: the population synthesis fits 
indicate that there is little dust attenuation in the op- 
tical, and at z = 4.42 the Chandra hard band probes 
very hard X-rays (11 — 44 keV) that are largely unaf- 
fected by dust in the line of sight. The integrated in- 
trinsic stellar luminosity, from the best-fit ESF model, 
is Lstars — 2.0 X IO^^Lq, which constrains the ratio of 
the luminosity of the AGN to the stellar component of 
HDF-G4 to Lagn I Lstars ^0.1; Cl is clearly starburst 
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Fig. 4. — Schematic of the merger simulation that best matches our observations. The top six panels show the stellar surface mass 
density looking down on the orbital plane. They are 60 h-^ kpc on a side , the time relative to the start of the simulation in /i ^ Gyr is 
given in the upper right hand corner, and the black curve shows the path of the center of the smaller component (C2). The bottom three 
panels show the remnant spheroid viewed from three different projections. 
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4. COMPARISON TO HYDRODYNAMICAL SIMULATION 

In order to gain some insight into the nature of the 
HDF-G4 system, we have employed hydrodynamic simu- 
lations and designed a model for the encounter. Our sim- 
ulations were preformed using an updated version of the 
publicly available N-Body/SP H (Smoothed P article Hy- 
drodynamics) code Gadget2 (|Springe]|l2005[ ). They in- 
clude star formation, supernova feedback, and black hole 
acc retion; for a de t ailed d escription of our methodology, 
see ISpringel et ahl (|2005a 'fR). Using this code, we track 
the interaction and merger of two disk gala x ies w hose 
properties are scaled, as in [Robertson et all ()2006[ ). to 
be appropriate for a redshift of z = 5. 

The progenitor galaxy models are motivated by the 
stellar mass and size of the observed components of the 
HDF-G4 system. In particular, the observed relative 
fluxes of CI and C2, assuming a constant MjL ratio, 
imply a 2:1 interaction. This is reproduced with progen- 
itor disks initialized with circular velocities of V200 — 320 
and 200 km s^^ respectively, where V200 is the Keplcrian 
circular velocity at a mean overdensity of 200pc- Both 
models have baryon fractions of 5% and initial gas frac- 
tions of /g = 0.8, which is consistent with gas fractions 
observed in z ^ 2 disks (Erb j3t al. 200 a), and is fur 



thcr motivated by models ( Hopkins et .al][2007c| ) to ex- 
plain large black hole to host stellar mass ratios at z ^ 2 
(jPeng et al.ll2006f ) . They are realized with 10^ and 5 x 10^ 



dark matter, and 1.2 x 10^ and 6 x lO"' baryonic particles 
respectively. 

Our high redshift merger model assumes that HDF-G4 
is witnessed soon after the first passage of the two galax- 
ies. This is based on the appreciable projected physical 
separation that implies a lower limit of 7 kpc on the in- 
ternuclear distance, regular appearance, and bridge of 
material that connects the two components. The bridge 
in particular, while potentially reflective of the irregular 
rest-frame UV m orphology of star-forming systems (e.g., 
iLaw et aLll2007[ ). is also a characteristic feature of galaxy 
interactions at low redshift and provides the strongest 
constraints on the mutual orbit. A prograde encounter 
for either galaxy can be ruled out based upon the ab- 
sence of the symmetric tidal features that these interac- 
tions produce. More likely is a fairly direct encounter in 
which the galaxies interpenetrate at first passage and ma- 
terial is strewn out to form a bridge between them. Af- 
ter several trials, the best fitting orbit was parabolic - in 
agreement with likely orbital geometries from cosraologi - 
cal simulations (|Bensonll2005HKhochfar fc Burkertll2006f) 
- with a small perigalacticon of 1 kpc and both disks ori- 
ented nearly perpendicular to the orbital plane. 

Snapshots of the stellar mass distribution, including 
three different projections of the remnant, at different 
stages of the encounter are shown in Figure [H W e treat 
the stellar particles as iBruzual fc CharlotI Hool) SSPs 
with ages and metallicities according to their formation 
time. Stars initialized with the progenitor disks are as- 
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Fig. 5. — The rest-frame 1500A surface brightness of the same 
snapshot in Figure |4] from the viewing angle that best matches the 
HDF-G4 system. This corresponds roughly to the observed i— band 
a.t z = 4.42. This image has the same scale and resolution as Fig- 
ure[T] and has been convolved with the ACS PSF. The black circle, 
as in Figure[T] represents l"or 7 kpc projected physical separation, 
and the red contours show the isophotal shapes. Note that the sim- 
ulations are able to roughly reproduce the rest— frame morphology 
of HDF— G4, with a bridge between the two components. 

sumed to have formed 0.7 Gyr before the start of the 
simulation with low metallicity Z = 1 x 10~^. We ne- 
glect absorption and scattering by dust because the stel- 
lar population synthesis analysis in § |3] indicates that 
HDF-G4 is virtually dust-free. The best match to the 
observed optical morphology (rest-frame UV A 1500 A) 
of HDF-G4 was att = 0.17 Gyr after the start of the sim- 
ulation, and is shown in Figure O The snapshot shows 
the rest-frame L ■> surface brightness at the resolu- 
tion of the observations in Figure [T] for z — 4.42, and 
convolved with the ACS PSF. The stellar mass at this 
snapshot is 2.9 x 10^° Mq with a total SFR of 59 Mq 
yr~^ ; both within the observational constraints given 
the uncertainty in the IMF and relative M/L ratios of 
CI and C2. The synthetic SEDs of the two components 
over the physical scales in the ACS image, shown in Fig- 
ure [51 are also a close match to the observations. Finally, 
the total luminosity from accretion onto the black hole is 
Lagn = 8 X 10® Lq, which is well within the upper limit 
of Lagn < 2 x IO^^Lq imposed by the observations. 
The ability of simulations to reproduce the observable 
features of this object supports our interpretation that 
HDF-G4 is an early-stage merger. 

Finally, we note as a caveat that the isophotal shapes in 
Figure [5] are somewhat more elongated than those in the 
optical imaging data. Though our population synthesis 
modeling suggests that the overall effects of dust in the 
HDF-G4 system are small, because this is the rest-frame 
UV, even small amounts of dust in these elongations 
could substantially change their appearance in simulated 
observations. The structure of these features is also sen- 
sitive to the specific orbital geometry and the detailed 
structure of the progenitor disks, both of which may not 
be precisely captured in our modeling. The bridge, on 
the other hand, is produced by ram-pressure stripping of 
gas when the progenitors interpenetrate. Furthermore, 
if there were small amounts of dust in the system, this 
stripping would not necessarily carry the dust along with 
it. As a result, though the detailed isophotal shapes of 
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Fig. 6. — Simulated fluxes at z = 4.42 versus our observations of 
the HDF— G4 system for both components over the physical scales 
in the ACS image. The squares and solid line are the observed and 
simulated SED for CI, and the triangles and dashed line are the 
observed and simulated SED for C2. 

the two components are sensitive to the initial condi- 
tions of our simulations, we believe that production of a 
UV-bright bridge between the two components is generic 
for any encounter with a small impact parameter, and is 
therefore the best evidence for an ongoing merger and 
the applicability of our modeling. 

5. DISCUSSION 

In Figure[71 we use our simulation to predict the future 
evolution of the IIDF-G4 system. This indicates that 
the two components will finally coalesce, with a peak in 
the starburst and AGN luminosity, around a redshift of 
z sa 3.5. At this time, the IIDF-G4 sy stem will be a 
Lboi « 4 X 10^^ Lq quasar. Based on the iHopkins et al.) 
(2006iQ) quasar luminosity function models, we would 
expect there to be of order a couple such systems in the 
GOODS-North field. Although we cannot robustly con- 
strain the space densities of these mergers, it is conser- 
vatively consistent with the expected counts. 

Our simulations also suggest that by a redshift of 
z « 2.5, IIDF-G4 will have a spheroidal morphology 
(see Figure 111), a stellar mass of M* = 7.4 x 10^° Mq, 
and a SFR of 8.8 Mq yr~^ ; a slowly star-forming 
but not entirely passive elliptical. Its colors will be 
red, with (J — K)ab ~ 1-2, which is close to the 
color cut o f (J — K)ar > 1-3 for Distant Red Galax- 
ies (PRCs: [van Dokkum et al]|2003f ). The re mnant also 
has a small effective radius - measured as in ICox et al] 
(j2006f ) - relative to local spheroids of the same mass, with 
Re = 1.4±0.17 kpc and an average inner stellar surface 
mass density within of CTg = (0.6 ± 0.03) x 10^° Mq 
kpc^^, where the errors take into account projection ef- 
fects. This makes it similar to the compa ct, quiescent 
high r edshif t spheroids a t z ~ 2 observed bv lLabbe et al.l 
(|2005[ ) and IZirm et al.l (|2007t) . who find effective radii 
of Re ~ 0.5 — 1.1 kpc, inner stellar surface mass den- 
sities of CTe ~ 0.7 — 4 Mq kpc~^, and stellar masses of 
- 5-9 X 10^° Mq. 

The existence of such an object further lends obser- 
vational support to galaxy formation models which ar- 
gue that mergers dominate the high- redshift buildup 
of spheroids and black holes (jHopkins et al.l l2007dl : 
ICroton et al.l l2006^. It is difficult, on the other hand, to 
reconcile this with models which argue that spheroid and 
black hole formation at high redshift is driven by direct 
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Fig. 7. — The future evolution of the SFR (top panel) and the 
AGN luminosity (top panel) of the HDF— G4 system as predicted 
by our simulations. This suggests that HDF— G4 will be a quasar 
at 2 Si 3.5, and a quiescent spheroid a,t z ^ 2.5. 



gas collapse /cooling (iGranato et al.ll20d^ or disk insta- 
bilities (e.g-. lBower et al.ll2006f l. In the latter model, for 
example, mergers contribute only 0.1% of bulge and 
BH growth at these redshifts, which would predict sys- 
tems like HDF-G4 should be too rare to observe in small 
fields, should be already "dry" (i.e. spheroid-dominated, 
with much lower current SFRs), and should have pre- 
existing large BHs from their disk phase (given the disk 
masses and gas fractions necessarily present at these red- 
shifts) which would, given gas is clearly still present 
(from the observed SFR), inescapably make the system a 
(~ 3 X 10^^ Lq) quasar at the observed time. Therefore, 
it will be interesting to see if future observations of wider 
fields reveal HDF-G4 to be anomalous, or representative 
of a typical, albeit brief, phase of LBG evolution. 

6. CONCLUSION 

We present the rest-frame UV through NIR SED 
for HDF-G4, an interacting Lyman Break galaxy in 



GOODS-North with a spectroscopical ly confirmed red- 
shift of z = 4.42 (Ste idel et al.lll999f ). using data ob- 
tained from AGS, WIRCam, and IRAC. The two objects 
in this system - HDF-G4 and its previously unidenti- 
fied companion - are both ^435 band dropouts, have 
similar Vgoe — 1775 and 1775 — zgso colors, and are sep- 
arated by 1", which at z = 4.42 corresponds to 7 kpc 
projected nuclear separation, and a bridge of material be- 
twe en them. We apply stella r population synthesis mod- 
els (iBruzual fc Char loll 120031 to the SED of HDF-G4, 
and find a best-fit population with solar metallically, an 
exponential star formation history, with e-folding time 
To = 440 My, stellar age = 720 Myr, stellar mass 
= 2.6 X lO^^Mg, and dust extinction E{B -V) = 
with a reduced x^/d.o.f = 1.46. The observed stellar 
population, combined with the SED, system morphol- 
ogy, and projected nuclear separation, are used to con- 
strain a model of the HDF-G4 system using a hydrody- 
namical simulation. This analysis suggests that HDF- 
G4 is the potential progenitor of a z 3.5 quasar with 
L PS 4 X 10^^ Lq, and a compact {Re = 1.4 kpc) qui- 
escent z 2.5 spheroid c onsist ent with the populatio n 
observed bv iLabU et"all (|2005l ) and IZirm etall (|2007f l. 
Furthermore, the existence of such an object supports 
galaxy formation models in which major mergers drive 
the high redshift buildup of spheroids and black holes 
(e.g.. lHopkins et al.ll2007dl ). 

We thank the referee for their helpful comments. This 
work is based on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, un- 
der NASA contact 1407. CFHT observations are sup- 
ported through the Taiwan CosPA project. 
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